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We present a comprehensive study of the discretized modes of an atomic gas in different conditions 
of confinement. Starting from the equations of hydrodynamics we derive a closed equation for the 
velocity field, depending on the adiabatic and isothermal compressibilities and applicable to different 
dimensions and quantum statistics. At zero temperature the equation reproduces the irrotational 
behavior of superfluid hydrodynamics. It is also applicable above the critical temperature in the 
collisional regime, where the appearence of rotational components in the velocity field is caused by 
the external potential. In the presence of harmonic trapping, a general class of analytic solutions is 
obtained for systems exhibiting a polytropic equation of state, characterized by a power law isoen- 
tropic dependence of the pressure on the density. Explicit results for the compressional modes are 
derived for both Bose and Fermi gases in the pancake, cigar as well as in the deep 2D and ID regimes. 
Our results agree with the analytical predictions available in the literature in some limiting cases. 
They are particularly relevant in ID configurations, where the study of the collective frequencies 
could provide a unique test of the achievement of the collisional regime at finite temperature. 

PACS numbers: PACS numbers 


I. INTRODUCTION 

In recent years the study of low-dimensional quantum 
gases has been the object of significant experimental and 
theoretical research [1, 2]. The main interest is due to 
the increase of quantum correlations and fluctuactions 
caused by the reduced dimensionality which make the 
properties of these systems significantly different from 
their 3D counterpart [3]. 

In two dimensions, thermal fluctuations rule out the 
existence of long range order at finite temperature 
(Hohenberg-Mcrmin-Wagner theorem) [4, 5] and novel 
aspects of superfluid phenomena, like the Berezinskii- 
Kosterlitz-Thouless transition, take place. In compress¬ 
ible ID systems the occurrence of long range order is 
ruled out by quantum fluctuations even at T = 0 [1, 6]. 
The goal of this paper is to focus on the collective os¬ 
cillations of a low dimensional atomic gas. While some 
of the results discussed in this paper have been already 
derived in the literature, a unified description covering 
both Bose and Fermi statistics at zero as well as at finite 
temperature in different regimes of axial and radial trap¬ 
ping is still missing. Our analysis can help in providing 
useful links and comparisons among different regimes of 
experimental and theoretical interest. 

Our investigation is based on the hydrodynamic for¬ 
malism applicable to various regimes, including the su¬ 
perfluid regime at zero temperature and the collisional 
regime in the non superfluid phase above the critical tem¬ 
perature. To this purpose, we will derive a general for¬ 
mulation of the hydrodynamic equations in presence of 
an external trapping potential in terms of the velocity 
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field. This approach explicitly points out the irrotational 
vs rotational nature of the solutions and the role of the 
equation of state. It reduces to a simplified form when 
the equation of state exhibits a polytropic dependence 
on the density, allowing for an important class of ana¬ 
lytic solutions. Our results are also compared with the 
predictions for the frequencies of the normal modes in 
the collisionless regime. 

The paper is organized as follows: 

In Sec. II we introduce the concept of polytropic equa¬ 
tion of state, both in terms of the pressure and, at T = 0, 
of the chemical potential. The values of the polytropic 
coefficient q , fixing the power law isoentropic dependence 
of the equation of state on the density, are discussed in 
some relevant cases of different dimensionality for both 
Bose and Fermi gases at zero as well as at finite temper¬ 
ature. We will consider geometrical configurations corre¬ 
sponding to pancake and cigar profiles, whose equation 
of state can be derived starting from the 3D equation of 
state, employing the local density approximation along 
the directions of the confinement. We will also discuss 
the case of deep 2D and ID configurations where the mo¬ 
tion of the gas is instead frozen along the directions of 
the confinement. 

In Sec. Ill we derive a general equation obeyed by the 
velocity field associated with the collective motion, start¬ 
ing from hydrodynamic theory in the presence of external 
trapping. This equation holds both for a superfluid gas 
at zero temperature, where the macroscopic dynamic be¬ 
havior is described by the hydrodynamic theory of super¬ 
fluids, and in the normal phase, in the collisional regime, 
characterized by the condition uit -C 1 , where r is a typ¬ 
ical collisional time. We will show that the equation for 
the velocity field acquires a particular simple form if the 
equation of state of the gas exhibits a polytropic depen¬ 
dence on the density. 

In Secs. IV, V, VI and VII we derive the discretized 
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oscillation frequencies in the presence of additional har¬ 
monic trapping in the plane (in the case of 2D configu¬ 
rations) and along the axis (in the case of ID configura¬ 
tions). Special emphasis is given to the lowest breathing 
modes whose frequency, differently from the divergency 
free solutions, is explicitly sensitive to the value of the 
polytropic coefficient and can exhibit a temperature de¬ 
pendence even in the hydrodynamic regime. 

II. POLYTROPIC EQUATION OF STATE 

A useful feature, which allows for a significant simpli¬ 
fication of the theoretical analysis, is the fact that in 
several configurations, corresponding to different dimen¬ 
sional, interaction and temperature regimes and quan¬ 
tum statistics, the equation of state of the gas exhibits 
a simple power law dependence on the density n in the 
form 

P(n, s) = n q p{s) , (1) 

where P is the pressure of the gas, q is called the poly¬ 
tropic coefficient, s is the entropy per particle and p(s) 
is a function of s, fixed by the thermodynamic behav¬ 
ior of the gas at zero (s = 0) as well as at finite tem¬ 
perature. As we will show, the polytropic coefficient q 
plays a crucial role in determining the frequency of the 
discretized collective oscillations in the presence of har¬ 
monic trapping. Viceversa, the value of the function p(s) 
determines the adiabatic sound velocity, according to the 
thermodynamic law 

me 2 = (§Q _ = ^ 9_ V(s) = ■ (2) 

Starting from Eq. (1) and using the thermodynamic 
relation ( dU/dn ) g = P/n 2 , one can write the following 
relationship between the energy per particle and the pres¬ 
sure of the gas: 


At T = 0 the scaling law (1) can be also written in the 
form 

T = 0) = = 0) (4) 

where p is the chemical potential of the gas related to the 
pressure by the Gibbs-Duhem thermodynamic relation 


When applied to configurations of lower dimensions the 
pressure and the density entering (1) should be replaced 
by the corresponding 2D or ID pressure and density of 
the gas. 

The parametrization (1) for the equation of state ap¬ 
plies to an important class of physical systems. The val¬ 
ues of the corresponding polytropic index q are explicitly 
derived in Appendix A and are reported in the Tables I 
and II. For example, in a 3D Fermi gas at unitarity [7] 
(as well as in the ideal 3D Fermi gas and in the ideal 3D 
Bose gas above the critical temperature) simple dimen¬ 
sionality arguments permit to identify the value q = 5/3 
for the polytropic coefficient. The 3D weakly interacting 
Bose gas [8], at zero temperature, is instead character¬ 
ized by the different value q = 2. On the other hand, Eq. 
(1) does not hold for the ideal Bose gas below the con¬ 
densation temperature, since this system has an infinite 
compressibility and its pressure vanishes in the limit of 
zero temperature. 

In the presence of axial or radial confinement one 
should distinguish between the case when the 3D equa¬ 
tion of state can be still applied locally within the local 
density approximation (LDA) and the case where the mo¬ 
tion is instead frozen along the directions of confinement. 

The first case includes the so called pancake geometry 
(where the equation of state can be expressed in terms 
of the 2D density and pressure) and the so called cigar 
geometry (where the equation of state can be expressed 
in terms of the ID density and pressure). Both configura¬ 
tions have been realized experimentally both for bosons 
[9, 10] as well as for fermions [11-13] and are character¬ 
ized by different values of the polytropic coefficient (see 
Table I). The second case includes the deep 2D and ID 
configurations where the thermodynamic behavior can¬ 
not be derived starting from the 3D equation of state 
in the LDA. These configurations are particularly inter¬ 
esting from the many-body and thermodynamic point 
of view. In 2D they exhibit the Berezinski-Kosterlitz- 
Thouless transition [14], while in the ID case for bosons, 
at zero temperature, they are described by Lieb-Lininger 
theory [15]. Some of these low dimensional regimes are 
also characterized by well defined values of the polytropic 
coefficient which are reported in Table II. 

It is finally interesting to compare the adiabatic sound 
velocity c^£> of 3D uniform systems with the sound veloc¬ 
ity c for the lower dimensional regimes reported in Table 
I. At T = 0, for equal values of the chemical potential 
(which, in the case of the LDA regimes, would correspond 
to equal values of the central density), one finds the fol¬ 
lowing relationship 


dP = ndp + sdT (5) 

where s = ns is the entropy density. The sound velocity 
at zero temperature, where dP = ndp , can be rewritten 
in terms of the chemical potential as 

me 2 = ={q-l)p . (6) 


2 i • V' 7 

c 3 D ?3D - 1 

where q^o is the polytropic coefficient of the equation 
of state of 3D uniform matter. In the special case of 
a cigar trap (see Table I), the above equation gives the 
value c 2 /c % D = 1/2 for a weakly interacting Bose gas and 
c 2 /c\ D = 3/5 for the unitary Fermi gas. The above reduc- 
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TABLE I: Polytropic index q for the weakly interacting Bose 
gas at T = 0 and for the unitary Fermi gas in 
different LDA regimes. The results for the unitary 
Fermi gas hold also for the ideal Fermi gas at all 
temperatures and for the ideal Bose gas above T c . 



Bose gas (T 

= 0) Unitary Fermi gas 

3D uniform 

2 

5/3 

pancake (LDA) 

5/3 

3/2 

cigar (LDA) 

3/2 

7/5 


TABLE II: Polytropic index q for a Bose gas for different 2D 
and ID regimes in the presence of tight 
confinement. The 2D mean field value q = 2 
holds also for the interacting Fermi gas in the 
BCS and BEC limits. 



T = 0 high T 

2D mean held 

2 2 

ID mean held 

2 3 

ID Tonks-Girardeau 

3 3 



tion factors of the sound velocity were first theoretically 
derived by Zaremba [16] and by Capuzzi et al. [17]. They 
were confirmed experimentally by Andrews et al [18] and 
Joseph et al. [19] for cigar Bose and Fermi gases, respec¬ 
tively. In the classical regime of high temperature, where 
P = nkT , one instead find the result 


where the sound velocities are calculated at the same 
temperature. 

III. HYDRODYNAMIC EQUATIONS IN THE 
PRESENCE OF EXTERNAL TRAPPING 


and Fermi superfluid gases, as well as to strongly inter¬ 
acting superfluids like 4 He [1], They also apply above 
the critical temperature Tc for superfluidity, in the colli- 
sional regime where lot <C 1, with r the collisional time, 
and the equation of state can still depend in a crucial way 
on quantum statistical effects [20]. For the ideal Fermi 
gas they were derived in the collisional limit in [21] at T = 
0. For finite temperature, below Tc, the hydrodynamic 
equations should be instead generalized to include the 
coupled description of the normal and superfluid compo¬ 
nents (Landau’s two-fluid hydrodynamic equations), al¬ 
lowing for the propagation of first and second sound (see, 
for example, [1]). In the presence of harmonic trapping 
the hydrodynamic equations have proven quite success¬ 
ful at zero temperature in predicting the frequencies of 
the collective oscillations [22] , in fair agreement with the 
experimental findings. 

The hydrodynamic equations include the equation of 
continuity 

+ V • [n 0 (r)v(r, t)] = 0, (9) 

the equation for the entropy density 

9S g’ ^ = -V [so(r)v(r, *)] , (10) 

and the Euler equation for the current density j(r,f) = 
mn 0 (r)v(r, t) 

= - [VP(r, t) + n{ r,f)V V ext {r)\ (11) 

In the above equations no and so are the equilibrium val¬ 
ues of the density and of the entropy density, respectively 
and, for simplicity, we have considered the limit of small 
amplitude oscillations and small velocities. At zero tem¬ 
perature, where the equation for the entropy identically 
vanishes, the hydrodynamic equations reduce to two cou¬ 
pled equations for the density and the velocity field. 

Starting from the above equations and using suitable 
thermodynamic relations, it is possible to derive the fol¬ 
lowing equation for the stationary solutions of the veloc¬ 
ity field characterized by the time dependence v(r, t ) = 
v(r)e _lwt (see Appendix B): 


In this Section we discuss the hydrodynamic (HD) 
equations for the velocity field describing the collective 
motion of the gas in the presence of an external confining 
potential and in the dissipationless (entropy conserving) 
regime. Even if these equations have a classical shape, 
quantum mechanical effects deeply affect their solutions, 
being responsible for important changes in the behavior 
of the equation of state and, in particular, in the value of 
the polytropic coefficients, as discussed in the previous 
Section. These equations describe correctly the dynamic 
behavior of an interacting system at zero temperature, 
where they coincide with the irrotational hydrodynamic 
equations of superfluids and apply to interacting Bose 



(V-v) 


+ ( 7 -l)(Vb elt )(V-v) 


+ V (v • W ext ) (12) 


where the adiabatic coefficient 



provides the ratio between the isothermal and the adi¬ 
abatic compressibilities. Equation (12) holds for con¬ 
figurations of different dimensionality and for arbitrary 
external potentials compatible with the applicability of 
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the local density approximation along the direction of 
the velocity flow. It shows that the eigenfrequencies u> 
of the collective oscillations can be determined once the 
adiabatic and the isothermal compressibilities are known. 
These quantities depend on the nature of the system, on 
temperature and, of course, on the dimensionality of the 
configuration. Result (12) shows explicitly the emergence 
of rotational components in the velocity field v caused by 
the presence of the external potential (term proportional 
to (7 — 1)). This term exactly vanishes if the adiabatic 
and isothermal compressibilities coincide (7 = 1 ), a con¬ 
dition ensured only at T = 0. In this limit Eq. (12) 
actually coincides with the hydrodynamic equation of ir- 
rotational superfluids. To our knowledge, the equations 
of linearized hydrodynamics, in the presence of an ex¬ 
ternal potential, have never been derived so far in the 
general form (12). Their derivation represents one of the 
main results of the present work. 

Starting from the equations of hydrodynamics it is also 
possible to derive (see Appendix B) the equation 


solutions that will be discussed in the following Sections 
where we will consider an external potential of the shape 

KxtOc, y) = (u 2 x x 2 + u 2 y y 2 ) (18) 

for the investigation of the collective oscillations of 2D 
configurations and of the form 

V ext (z) = muj 2 z 2 (19) 

in the case of ID configurations. 

IV. SURFACE MODES 

For divergency free oscillations (V • v = 0), also called 
surface modes, the hydrodynamic equation ( 12 ) takes the 
simplified form 


dT(r,t) T (dP\ V • v(r,t) 

dt c v \dT) n n 0 [ ’ 

for the time dependence of the temperature of the gas, 
where c v is the specific heat at constant volume. Eq. 
(14) shows that, for divergency free solutions (V • v = 0), 
like the surface modes (Sec. IV), the temperature of the 
trapped gas is constant during the oscillation. 

For systems obeying the polytropic behavior (1) the 
equation for the temperature takes the simplified form 

^M=-(?-l)TVw(r ,t) (15) 

which follows directly from Eq. (3) and from the ther¬ 
modynamic relation c v = (dU/dT) n . Equation (15) co¬ 
incides with the result derived in [20] for a 3D ideal Bose 
gas at T > T c , where q = 5/3. 

For systems obeying the poly tropic law (1), the equa¬ 
tion for the velocity field ( 12 ) takes the further simplified 
expression (see Appendix B) 

mw 2 v = - (jjPj V (V • v) + (q - 1) (Wezt)V • v 

+ V(vVU (16) 

in terms of the adiabatic compressibility and the poly¬ 
tropic coefficient q. Equation (16) coincides with the re¬ 
sult derived in [20] for a 3D ideal Bose gas for T > T c 
with q = 5/3. At T = 0, where 7 = 1, it reduces to 

mcu 2 v = -(q- 1)V [/z(n)V • v] + V [v • VV ext ] (17) 

where we have used the T = 0 expression (4) for the 
chemical potential y(n). 

In the form (16) (and (17) at T = 0) the equation 
for the velocity field admits a useful class of analytical 


w 2 v = V (v • VV ext ) (20) 


independent of the equation of state, quantum statistics 
and temperature. 

In the case of 2D isotropic trapping ( uj x = oj y = w_l) 
an important class of solutions of Eq. (20) is given by 
the irrotational choice [ 22 ]: 


v(rj_,0) oc V 


\ m \ ±im<t> 

/ e , 


( 21 ) 


where m is z-th component of angular momentum. The 
resulting eigenfrequencies take the form 

w 2 (m) = uj 2 1 _\m\ . ( 22 ) 

Important examples of divergency free oscillations, sys¬ 
tematically investigated in harmonically trapped atomic 
gases, are the dipole (ojd = wi) and the quadrupole 
(ojq = -\/ 2 w_l) oscillations. The former case is model in¬ 
dependent, due to Kohn’s theorem. The latter case is 
particularly interesting because the predicted frequency 
differs from the collisionless ideal gas value 2u)± and its 
investigation can then be used as a direct test of the 
achievement of the hydrodynamic nature of the regime 
considered. 

In the case of anisotropic potentials ( lv x ^ uj v ) an 
important example of divergency free oscillations is de¬ 
scribed by the ansatz: 

v(x,y) ocV(xy) (23) 

yielding the result 

4 =^+^ ( 24 ) 

for the collective frequency. This is the so-called scis¬ 
sors mode corresponding to an oscillating rotation of the 
atomic cloud in the x-y plane. It was first predicted by 
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Guery-Odelin and Stringari (1999) [23] and experimen¬ 
tally observed in both 3D Bose [24] and Fermi superflu¬ 
ids [25] . This mode preserves the shape of the gas during 
the oscillation and for this reason is independent of the 
equation of state. It has the same value in 3D as well 
as in 2D configurations described by the hydrodynamic 
formalism. 


V. COMPRESSIONAL OSCILLATIONS IN 
PANCAKE AND 2D ISOTROPIC TRAPS 

Let us first consider the superfluid zero temperature 
regime for which, as already pointed out in the Section 
III, the irrotationality condition for the velocity field 
holds. For systems whose equation of state is character¬ 
ized by the polytropic density dependence (4), the veloc¬ 
ity field obeys the hydrodynamic Eq. (17). For isotropic 
2D traps one can easily check that this equation is solved 
by the following irrotational ansatz for the velocity field 
[26]: 


v(rj_, </>) oc V 


{r 2 ? + 


(25) 


where n fixes the number of the radial nodes of the den¬ 
sity modulations occurring during the oscillation and to 
is z-th component of angular momentum. The resulting 
hydrodynamic eigenfrequencies take the form 

uj 2 (n, to, q) = w5_[2n + |to| + 2n(q — l)(n + |m|)]. (26) 

Result (26) can be applied to a rich variety of physical 
situations whose equation of state is simply incapsulated 
in the polytropic q index and it differs from the predicted 
values of the ideal gas model in the collisionless regime 

u) c i(n, m) = wj_(2n + |to|). (27) 

Moreover, result (26) should be compared with the hy¬ 
drodynamic dispersion law holding in the presence of a 
3D isotropic harmonic potential (lo x = lo v = ui z = Uho ) 

w 2 (n, l, q) = w 2 D [2n + l + (q — l)(2n 2 + 2 nl + n)] (28) 

where l is the angular momentum. Eq. (28) reduces 
to the results of [22] and [21, 27, 28] in the case of the 
weakly interacting Bose (q = 2) and Fermi (q = 5/3) 
gases, respectively. Result (26) reduces to the dispersion 
(22) of surface modes in the case n = 0. 

The most important solution accounted for by (26) is 
the 2D breathing or monopole mode (n = 1 and to = 0), 
for which one finds the result 


WAf(g) = 2qu\. (29) 


frequency (29) of the 2D breathing mode differs from the 
one in the 3D isotropic case (see Eq.(28), n = 1 and l = 0) 
which yields the value Loj^(q) = (3 q — l)w 2 0 . 


For the unitary Fermi gas in the pancake regime (q = 
3/2), one gets the value lom = V3oj±. Actually this 
value turns out to be independent of temperature, as a 
consequence of the fact that, at unitarity, the scaling so¬ 
lution for the breathing mode is exactly satisfied by the 
hydrodynamic equations at all temperatures [29]. No¬ 
tice that this result differs from the value lom = 2 coho 
holding for the unitary Fermi gas with 3D isotropic trap¬ 
ping (q = 5/3). In this latter case an exact solution of 
the time dependent Schrodinger equation is available for 
the monopole breathing mode, independent not only of 
temperature but also of the amplitude of the oscillation 
and holding in all collisional regimes (scale invariance) 
[30]. A similar situation holds also for the 2D regime 
of a weakly interacting Bose gas (q = 2) in the pres¬ 
ence of isotropic trapping where one also finds the result 
lom = 2w_l, independent of temperature and of the ampli¬ 
tude of oscillations [31, 32]. In the case of the 2D Fermi 
gas the result lom = 2u;j_, following from the q = 2 value 
of the polytropic coefficient, holds both in the BEC and 
BCS regime and, with high accuracy, along the whole 
crossover, revealing an apparent scale invariance [33, 34], 
as proven experimentally in [35]. 


The temperature dependence of the frequency of the 
breathing mode takes instead place for the pancake 
weakly interacting Bose gas. At T = 0 (q = 5/3) one 
finds ojm = v/10/3w_l [22]. This result differs from the 
value obtained at high temperature, where the thermo¬ 
dynamic behaviour of the gas can be approximated by 
the ideal Bose gas {q = 3/2) and the hydrodynamic fre¬ 
quency takes the value % = V3ui±. The frequency of 
the monopole mode of the pancake Bose gas is then ex¬ 
pected to exhibit a temperature dependence. The tem¬ 
perature dependence of this mode was pointed out in the 
first experiments on the collective oscillations carried out 
at Jila [36] where, however, due to the small number of 
atoms, the system enters soon the collisionless regime for 
T > T c , rather than the collisional one. Achieving the hy¬ 
drodynamic condition lot <C 1 is in general a difficult task 
in weakly interacting Bose gases, due the small available 
values of the density of the gas. A simple estimate, based 
on the classical evaluation of the collisional time, yields, 
for frequencies of the order of the trapping frequency Uhot 
the condition [8] l m ph Rt, where l m ph = (tot) -1 is the 
mean free path. Here n is the local 3D de nsity, a = 87t a 2 
is the s-wave cross section, while Rt = \j2ksT I mLo\ 0 is 
the thermal radius. 


Differently from the surface modes discussed in the previ¬ 
ous Section, the frequency of the breathing mode depends 
explicitly on the polytropic index q and is consequently 
sensitive to the equation of state obeyed by the gas. The 


A summary of the frequencies of the breathing mode 
in the most relevant cases discussed in this section is 
presented in Tables III and IV. 
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TABLE III: Hydrodynamic frequencies of the breathing 
mode in different LDA regimes. The values 
reported for the Unitary Fermi gas hold also for 



the ideal Bose 

gas above T c . 


Bose gas (T = 

0) Unitary Fermi gas 

3D isotropic 

V^UJho 


pancake (LDA) 

y/W/&>± 

V3ui± 

cigar (LDA) 

y/5/2u z 

yi2/5u>z 


TABLE IV: Hydrodynamic frequencies of the breathing 

mode for a Bose gas in low dimensions and tight 
trapping regimes. The 2D mean field value 
lom = 2cuj_ holds also for the interacting Fermi 
gas in the BCS and BEC limits. 



T = 0 

high T 

2D mean field 

2u>± 

2u>± 

ID mean field 


2cj z 

ID Tonks-Girardeau 

to 

£ 

2io z 



VI. COUPLED COMPRESSIONAL AND 

SURFACE MODES IN ANISOTROPIC TRAPS 

In the presence of anisotropic trapping (ui x ^ u y ), the 
breathing and the quadrupolc oscillations discussed in 
the previous Sections are coupled. The coupling is ac¬ 
counted for by the ansatz 

v(x,y) = V(ax 2 +/3y 2 ), (30) 

where the relative value of the parameters a and /? have 
to be determined by solving the hydrodynamic equation 
(16). Since V (V • v) =0 the frequency of the coupled 
modes, for a given value of the polytropic coefficient q, 
are given by 

w) cy v^x ^y) 

J(q+ 1 ) 2 (w2 + W 2)2 _ 16 

^--■ ( 31 ) 

In the isotropic limit (ui x = u> y = wj_) Eq. (31) re¬ 
produces the results (29) and (22) for the monopole and 
the quadrupole, respectively, derived in the previous Sec¬ 
tions. In the case of the 2D mean field regime, where the 
polytropic coefficient takes the value q = 2, result (31) 
was derived in [37] and in [38] for bosons and fermions, 
respectively. 

Result (31) gives general predictions for the collective 
frequencies in 2D configurations for arbitrary values of 


the deformation of the trap. It is interesting to consider 
the limit of highly deformed 2D trapping potentials. For 
oj x *C Uy the lower solution takes the form 


JL = 


4 9 


(9 + 1 ) 


while the upper solution is given by 

= (9 + i) w y 


(32) 


(33) 


Of course a symmetric result takes place in the opposite 

limit oj x uiy. 

In the 2D mean field case (q = 2) (32) provides the 
result = -y/8/3u; x , in agreement with the findings of 
[37]. In the pancake regime (32) instead reduces to the 
frequency of the axial breathing mode of a cigar configu¬ 
ration (see next Section). 

Equation (33) instead coincides with the solution of the 
3D hydrodynamic equations for triaxial harmonic trap¬ 
ping [1] in the intermediate regime 


VII. COLLECTIVE FREQUENCIES IN CIGAR 
AND ID TRAPS 

In this Section we complete the discussion on the col¬ 
lective frequencies of harmonically trapped gases consid¬ 
ering one dimensional configurations. 

Similarly to the 2D case, two different regimes can be 
considered also in ID. The first one, called cigar regime, 
corresponds to systems described locally by the 3D equa¬ 
tion of state, but the radial confinement is enough tight to 
ensure a ID nature to the low energy dynamic behavior. 
This configuration is particularly suited to investigate the 
propagation of sound [18, 19, 39]. 

A second case is the deep ID (hereafter simply called 
ID) regime where the radial motion is frozen (for a gen¬ 
eral review on ID systems see, for example, [40, 41]). 
The role of correlations is particularly important in this 
regime and, in the case of Bose gases at T = 0, is de¬ 
scribed by Lieb-Lininger theory [15]. A key question in 
ID is the role of thermalization. From this point of view 
the comparison between experiments and the predictions 
of hydrodynamic theory (which assumes local thermal¬ 
ization at finite temperature) would be particularly in¬ 
sightful. Another important feature of ID systems is the 
absence of phase transitions at finite temperature [42] . 

The equations of hydrodynamics derived in Sect. Ill 
become particularly simple in ID, the velocity field being 
a function of the variable z and all the gradients acting 
only on the ^-direction. Looking for solutions of the form 

v = z k + a k - 2 Z k ~ 2 + ... (34) 

one finds that the hydrodynamic equation (16), in the 
presence of the poly tropic equation of state (1), admits 
simple analytic solutions in the presence of the harmonic 
trapping (19). 






















7 


At T = 0, where one can conveniently use Eq. (17) 
with the chemical potential given by (Al), the dispersion 
law takes the form 


w 2 (fc,g) = (k + l) 






(35) 


Equation (35) was first derived in [43] in the context of 
ID Bose gases. 

In the classical limit of high temperature, where one 
can use the equation of state Pm = ninksT^ the hydro- 
dynamic equation (16) gives rise instead to the following 
dispersion relation 

w 2 (fc,g) = (qk + 1 )uj 2 z (36) 

Result (35) should be also compared with the disper¬ 
sion 


w c z(/c) = (k + l)u z (37) 

holding in the collisionless regime of a non interacting ID 
gas. 

Let us first consider the case of the cigar unitary Fermi 
gas. In this case the polytropic coefficient is equal to 
q = 7/5 at all temperatures and the k = 1 hydrody¬ 
namic solution, corresponding to the lowest breathing 
(LB) mode, is equal to ujlb = -\JV2/hu z [44, 45] at both 
zero and finite temperature. The situation is different for 
the higher nodal modes {k > 1 ) where the hydrodynamic 
frequencies at T = 0 and in the classical limit are dif¬ 
ferent revealing an interesting temperature dependence 
[46] that was investigated experimentally in [47] in good 
agreement with the predictions of theory. The case of the 
cigar Bose gas is different. In fact in this case the poly¬ 
tropic coefficient depends on temperature, being equal to 
<7 = 3/2 at T = 0 and to q = 7/5 in the classical limit. 
As a result, the frequency of the k = 1 lowest breathing 
mode takes the value ujlb = -\/5/ 2uj z [22, 26] while in 
the high temperature classical limit one finds the smaller 
hydrodynamic value to lb = \fV2/5u> z . In the collision¬ 
less limit one instead finds the value 2 u> z . The T = 0 
value yj5/2w z was measured experimentally in [48] in 
excellent agreement with the prediction of theory. At fi¬ 
nite temperature these authors found that the frequency 
slightly drops below the low temperature limit. At even 
higher temperatures they observed a significant increase 
of ui, likely due to the breakdown of the hydrodynamic 
condition lot -C 1 . 

The comparison between the lowest compressional 
mode in the cigar geometry [k = 1) and the lowest so¬ 
lution (32) holding in the pancake geometry with highly 
anisotropic 2D trapping, allows for the non trivial rela¬ 
tionship 

qa a ar + l= Aqpancak * (38) 

Qpancake i 1 

between the polytropic coefficients in the cigar and pan¬ 
cake geometries. The relationship is confirmed by the 
results reported in Table I. 


The behavior of the collective frequencies in the deep 
ID regime is also very interesting. In the case of bosons, 
one should distinguish between the mean field case, where 
q = 2 at T = 0, and the Tonks-Girardeau (TG) limit [49], 
where q = 3. In the former case the frequency of the k = 1 
breathing mode takes the value lolb = V^3w z , while in 
the TG limit one finds ujlb = 2 io z . The behaviour of the 
lowest breathing frequencies, in the intermediate regimes 
between the mean-field and the TG limits, was theoret¬ 
ically investigated in [43] and experimentally observed 
in [50]. At high temperatures, both regimes predict the 
frequency ujlb = 2 lo z for the same k = 1 mode. One 
then concludes that, differently from the ID mean field, 
the frequency of the breathing mode in the TG regime 
is temperature independent. The result is not surprising 
since in the TG limit bosons behave like ID non inter¬ 
acting fermions where scaling invariance applies. It is 
worth noting, however, that the hydrodynamic frequen¬ 
cies of the higher nodal modes depend on temperature 
even in the TG limit. For example, the k = 2 mode has 
frequency uok -2 = 3 uo z at T = 0 and L 0 k =2 = V7lo z in 
the classical limit. The temperature dependence of the 
breathing mode in the ID mean field has been recently 
measured in [51]. 


VIII. CONCLUSIONS 

In this paper we have derived a unified description of 
the discretized collective oscillations of quantum gases in 
different conditions of trapping and dimensionality. 

A major result of this work is given by the derivation 
of a general hydrodynamic equation for the velocity field 
(Eq. (12)) depending solely on the adiabatic and isother¬ 
mal compressibilities of the gas. This equation takes a 
particularly simple form in the case of systems exhibiting 
a polytropic equation of state, characterized by a power 
law dependence of the pressure on the density (P oc n q ), 
for a fixed value of the entropy. The polytropic equation 
of state characterizes a significant class of many body 
configurations of either bosonic and fermionic nature. 

Analytic results for the collective frequencies have been 
calculated for pancake and cigar configurations, where 
the equation of state can be obtained from the 3D equa¬ 
tion of state using the local density approximation, as 
well as in the deep 2D and ID regimes where the motion 
is instead frozen along the direction of the confinement. 
Special emphasis is given to the lowest breathing mode, 
whose frequency has been shown to depend explicitly on 
the value of the polytropic coefficient < 7 , and to its cou¬ 
pling with the quadrupole oscillations in the presence of 
anisotropic external potentials. 

We have emphasized, in particular, the comparison be¬ 
tween the T = 0 results, characterizing the behavior of 
superfluids, and the high temperature behavior, both in 
the hydrodynamic and collisionless regime. In this re¬ 
spect an interesting perspective of the present work is 
given by the study of the collective frequencies in ID con- 







figurations at finite temperature which could provide a 
useful test of the equation of state of strongly interacting 
gases [52-54] as well as the achievement of the collisional 
(hydrodynamic) condition at finite temperature. 

Our results can be also used as a starting point to cal¬ 
culate the collective frequencies in systems whose equa¬ 
tion of state deviates from the polytropic law and the 
solution of the hydrodynamic equations can be handled 
using a perturbative approach. This procedure, first em¬ 
ployed to estimate the frequency shifts of dilute Bose 
gases caused by beyond mean field effects [55], was de¬ 
veloped in a systematic way by Astrakharchik [56] in a 
variety of configurations of different dimensionality and 
quantum statistics. It was recently applied by Merloti et 
al. [57] to explore the breakdown of scale invariance in 
a quasi-two-dimensional Bose gas due to the presence of 
the third dimension. 


Appendix A: Equation of state and polytropic 
coefficient 

In this Appendix we identify the polytropic coefficient 
q characterizing the equation of state (1) of a uniform 
gas in different low dimensional configurations and for 
different quantum statistics (see Tables I and II). 

Let us first consider the ideal situation of zero tempera¬ 
ture. As already discussed in Section II, a first interesting 
two dimensional regime is the so called pancake where the 
system keeps, locally, its 3D nature in the sense that the 
chemical potential is much larger than axial oscillator en¬ 
ergy y /§> hui z . In this case we can apply the local density 
approximation (LDA) along the axial direction and write 
the chemical potential in the form 

d o = M») + f (Al) 

which allows us to determine the ^-dependence of the den¬ 
sity profile. The value of y q coincides with the chemical 
potential y(n) calculated at z = 0 and is fixed by the nor¬ 
malization condition n 2 D = f dzn(r). The quantity yo 
plays consequently the role of the 2D chemical potential 
and exhibits an explicit dependence on the 2D density 
n, 2 D- This dependence characterizes the 2D equation of 
state and, for simplicity, we will omit the suffix 0 in y Q 
and use the simpler notation y(n 2 D ) to characterize the 
equation of state of the two-dimensional gas. 

The 2D equation of state of a pancake gas is easily 
derived in the case of dilute Bose gas where the 3D equa¬ 
tion of state has the form y = gn and the use of the LDA 
procedure (Al) gives the result [1] 

^M = (5^^£) 2/S . (A2) 

characterized by the value q = 5/3 for the polytropic 
coefficient. 

An analogous calculation can be carried out for a pan¬ 
cake Fermi gas at unitarity, where the 3D equation of 


state has the form y = £.B^-(37r 2 n) 2 / 3 , with £b the 
Bertsch parameter, yielding [1] 


3 /j3 

{n2D) = ( 27 t £§ w 2 — U2D 

' m 


1/2 


(A3) 


and hence the value q = 3/2 for the polytropic coefficient. 

Similar results can be obtained in cigar ID-like config¬ 
urations where the gas is harmonically trapped along the 
x — y directions. In this case, after integrating the radial 
3D Thomas-Fermi profile, one finds [1] 

y B {.n 1D ) = 2Huj±(ani D ) 1/2 (A4) 

and 

F, , f^Bh 2 \ 3/5 [15 2 \ 2/5 

d {ni D )= |-^-j \—nmu} ± n 1D J (A5) 


for the chemical potential of the cigar Bose and unitary 
Fermi gas respectively where nm = / dxdyn( r). From 
the above equations one derives the values q = 3/2 and 
q = 7/5 for the polytropic coefficients in the Bose and 
unitary Fermi cigars, respectively which are reported in 
the Table I. 

The results q = 3/2 and q = 7/5 for the polytropic co¬ 
efficients of the pancake and cigar Fermi gas at unitarity 
holds also at finite temperature. In fact the 3D result for 
the equation of state, given by Eq. (1) with q = 5/3, can 
be usefully rewritten in the form [46, 58] 

P 3D {x,T) = f p (x)-^- (A6) 

/\rjn 

and 

„ ( m rji\ _ fp( X ) _ fn(x) f A *7\ 

n 3 D\X,T) 3 — 3 (A7) 


where the temperature dependence of the 3D pressure 
and density follows from dimensionality arguments. We 
have here introduced the thermal wavelength A t = 
■\/2irh 2 /rnksT, the dimensionless functions f p (x) and its 
derivative f n {x) = f' p (x). The dimensionless parameter 
x = y/ksT, fixed by the ratio between the chemical po¬ 
tential and the temperature of the gas, determines the 
entropy per particle s according to the relationship 


s(x) _ 5 f p (x) 
k B 2 f n (x) 


(A8) 


The above results reflect the universality of the 3D uni¬ 
form Fermi gas at unitarity. 

By applying the local density approximation (Al) 
along the z-direction in the pancake geometry and along 
the x — y plane in the cigar geometry, after integration 
of the density and of the pressure profile along the same 
directions, from Eq. (A6) and Eq. (A7) one easily finds 
a new temperature dependence of the 2D and ID pres¬ 
sure and density for a given value of the entropy, fixed 
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by the ratio x = g,/k B T. For the pancake configuration 
one finds 


P2D(X2D,T) = -p-</ B 2 Fp D (x2D) (A9) 
y Tncu z 

and 

W 2 fl(^ 2 Di T) = TT\/ B 2 Pn D { X 2D) (A10) 

y 7TXUJ z 

while for the cigar ID configuration the result is [46] 

Pid(xid,T) = - -^-Fp D (xiD) (All) 

Ay TTWjJ | 

and 

tud(xid, T) = —3 - ^-F^ d (xid) (A12) 

Ay TTWjJ j 

holding at all temperatures. From the above results 
one easily finds, using Eq. (1) the values q = 3/2 and 
q = 7/5 for the polytropic coefficient in the 2D and 
ID cases, independent of temperature. In Eq. (A9) 
and (A10), we have defined the integrated functions 

F^(x 2 d) = dz'f n , p (x') with x 2 d = Ho/k B T, while 
in Eq. (All) and (A12) the analogous ID integrated 
functions are F/®(x\ D ) = 2n f^°° dr' ± r' ± f niP (x') with 
Xi d = Mo /k B T. The same results for the values of q hold 
also for the ideal Fermi gas as well as for the ideal Bose 
gas above T c where the thermodynamic functions can be 
written in the same form as for the unitary Fermi gas, 
the functions /„ and f p being of course different. 

So far we have considered the pancake and cigar ge¬ 
ometries where the local density approximation allows 
us to safely use the 3D equation of state locally. The 
situation is different in the opposite regime of tight axial 
or radial confinement where the motion is frozen to the 
lowest harmonic oscillator wave function along the tight 
directions. At zero temperature the condition for being 
in these deep 2D and ID regimes is given by /i -C Hoj z 
and n -C hui±, respectively. At high temperature the 
conditions are, instead, k B T <C fkj z and k B T -C fko±, 
respectively. At T = 0 these low dimensional regimes are 
easily described in the case of a weakly interacting Bose 
gas where the use of Gross-Pitaevskii theory yields the 
simple results for the 2D and ID equations of state 

M S ( n 2 £)) = g2Dn2D (A13) 


and 


{ni D ) = gioniD (A14) 


apart from unimportant additional constant terms. In 
Eq. (A13) the 2D coupling constant is given by [1] 


92D 



m a z 


(A15) 


and is fixed by the ratio a/a z between the 3D s-wave 
scattering length and the axial harmonic oscillator length 
a z = yjh/muj z , while in the ID case one finds the result 
[ 1 ] 


giD = 2ahuj± (A16) 

Result (A13) for the equation of state of the 2D Bose 
gas holds only in the limit of weak coupling (a -C a z ), 
ensuring the absence of quantum anomaly effects [59]. 
The ID equation of state (A14) instead holds provided 
niDd±/a 1. In both cases (also called 2D and ID 
mean field regimes) one identifies the value q = 2 for the 
polytropic coefficient. Both the 2D and ID mean field 
regimes have been achieved experimentally in [14] and 
[60, 61], respectively. 

In the ID case, extensive experimental and theoretical 
work has been done also to explore regimes beyond the 
mean field condition ni B a\/a 1 where the many-body 
properties of the gas are described by Lieb-Lininger the¬ 
ory. In the limit n\DO.\/a -C 1, the gas enters the limit of 
impenetrable bosons also called Tonks-Girardeau regime 
where the gas acquires a Fermi like behavior and its equa¬ 
tion of state exhibits a quadratic density dependence [49] : 

lx TG (n 1D ) = 71-2 (A17) 

In this regime the value of the polytropic coefficient is 

q = 3. 

In the case of interacting Fermi gases, the theoretical 
description of the lower dimensional regimes is more dif¬ 
ficult with respect to the Bose case. In two dimensions 
one can still identify a BCS-BEC crossover like in the 
3D case. However, in the two-body problem the pres¬ 
ence of a resonance always corresponds to the occurrence 
of a bound state, differently from what happens in 3D 
where this is ensured only for positive values of the 3D 
scattering length. Far from the resonance the chemical 
potential of the 2D Fermi gas is described by a linear 
dependence on the density, both on the BEC side, where 
the gas is described by a 2D system of bosonic molecules, 
and on the deep BCS limit, where the gas approaches 
the ideal Fermi gas behaviour. In both cases the value of 
the polytropic coefficient is q = 2. The same value of q 
characterizes the equation of state of the classical gas at 
high temperature. 

Tables I and II summarize the main results of this Sec¬ 
tion, reporting the values of the polytropic coefficients q 
in a class of interesting configurations. 


Appendix B: Hydrodynamic equations in the 
presence of an external potential 

In this Appendix we will derive the most relevant hy¬ 
drodynamic equations used in the paper to calculate 
the collective frequencies in the presence of external har¬ 
monic trapping. 
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In order to derive Eq. (12), we first take the time 
derivative of Eq. (11) yielding, in the limit of small am¬ 
plitude oscillations, 

”*"» w v =- v §; p - (!") ■ <Bi) 


The time derivative of the pressure can be written as: 


dP _ (dP\ dn fdP\ dT 
dt \ dn J T dt y dT ) n dt 


(B2) 


and requires the knowledge of the time derivative of the 
temperature. 

In order to calculate dT/dt it is convenient to rewrite 
the equation for the entropy density ( 10 ) in terms of the 
entropy per particle s = s/n. By using the equation of 
continuity (9), it is immediate to find the equation 


<9s(r, t) 
dt 


— v ( r , i)Vso(r) 


1 



v • Vn o 

n 


(B3) 


where, in deriving the second equality, we have used the 
thermodynamical relation 


T (dP\ 
Tds = c v dT — — 

applied to equilibrium ( dT = 0), where 


-^{&T) n \dr) s 


(B4) 


(B5) 


is the specific heat at constant volume. 

On the other hand, by considering s as a function of 
density and temperature, one can also write 


ds / ds \ dn / ds \ dT 
dt \ dn ) T dt \ dT) n dt 


(B 6 ) 


Using the equation of continuity (9) and thermodynamic 
relation (B3), one finally obtains the useful equation 


dT(r,t) 

dt 


T_ (dP\ V ■ v(r,t) 

c v \dTj n n 0 


(B7) 


for the time derivative of the temperature, holding also 
in the presence of external trapping. 


Using the equilibrium condition in the Eq. (11) 

VP ° = (S) Vn ° = (B8) 

and the thermodynamic relation 



relating the adiabatic and isothermal compressibilities, it 
is easy to recast the hydrodynamic equation (Bl) in the 
useful form 



(V-v) 


+ (7 — l)(V14it)(V • v) 


+ V (v • S7V ext ) (BIO) 


where 7 is the adiabatic coefficient defined by (13) and 
we have considered velocity fields oscillating as v(r, t) = 
v(r)e -iajt . 

Starting from Eq. (BIO), one can derive Eq. (16) in 
the case of a polytropic equation of state for which Eq. 
(2) holds. By considering, moreover, Eq. (5) at constant 
temperature, Eq. (Al) and Eq. (13), one can finally 
rewrite: 




V (V ■ v) 


+ q (V ■ v) VUearf 



(Bll) 


which plugged into Eq. (BIO) yields Eq. (16). 
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